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Abstract 

We formulate the massive domain wall fermions on anisotropic lattices. For the 
massive domain wall fermion, we find that the dispersion relation assumes the usual 
form in the low momentum region when the bare parameters are properly tuned. 
The quark self-energy and the quark field renormalization constants are calculated 
to one-loop in bare lattice perturbation theory. For light domain wall fermions, we 
verified that the chiral mode is stable against quantum fluctuations on anisotropic 
lattices. This calculation serves as a guidance for the tuning of the parameters in 
the quark action in future numerical simulations. 
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1 Introduction 



In recent years, considerable progress has been made in understanding chi- 
ral symmetry on the lattice. Domain wall fermions [1,2] and the overlap 
fermions [3,4,5,6] have emerged as two new candidates in the formulation 
of lattice fermions which have much better chiral properties than the conven- 
tional lattice Wilson or staggered fermions. Since chiral symmetry is so crucial 
to the theory of QCD, it is therefore desirable to use these new fermions if 
possible. However, due to their more expensive computational cost, this task 
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has not been fully accomplished, particularly for the study of multi-hadron 
states and hadrons with massive quarks. 

On the other hand, anisotropic lattices have been used to study heavy hadronic 
states and they proved to be extremely helpful in various applications. These 
include: glueball spectrum calculations [7,8], charmonium spectrum calcu- 
lations [9,10], charmed meson and charmed baryon calculations [11,12] and 
hadron-hadron scattering calculations [13,14,15,16]. Many of these studies in 
fact involve light quarks and chiral symmetry plays an essential role. It is there- 
fore more appropriate to study the light quarks using lattice fermions with 
better chiral properties. Using domain wall fermions with the physical four- 
dimensional lattice being isotropic, the heavy-light systems [17] and hadron- 
hadron scattering [18] have been studied. As we pointed out, anisotropic lat- 
tices can be very helpful in these studies. It is therefore desirable to study do- 
main wall fermions with the four-dimensional physical lattice being anisotropic 
and this will be the major purpose of this paper. Major possible applications 
of the anisotropic domain wall fermion approach that we have in mind are the 
heavy-light hadrons, exotic hadrons with light quarks and hadron-hadron scat- 
tering processes where the objects being studied are heavy and chiral property 
is crucial. 

In this paper, we study an explicit formulation of domain wall fermions on 
anisotropic four-dimensional lattices. We adopt the domain wall fermion action 
which is of Shamir type [2] . For the gauge action, we use the tadpole improved 
gauge actions [19,7] which have been used in glueball calculations. Similar to 
the case of Wilson fermions on anisotropic lattices, domain wall fermion action 
on anisotropic lattices contains more parameters than its isotropic counter- 
parts. Some of these parameters have to be tuned properly in order to exhibit 
a correct continuum limit. This problem is first analyzed in the case of free 
domain wall fermions on anisotropic lattices. We find that, in order to re- 
store the normal relativistic dispersion relation for the quark, parameters of 
the fermion action has to be tuned accordingly. Then, we compute the quark 
propagator in lattice perturbation theory to one-loop. Quark field and quark 
mass renormalization constants are obtained for various values of the bare pa- 
rameters. We also discuss the choice of the parameters which maintains good 
chiral symmetry. This perturbative calculation serves as a guidance for further 
non-perturbative Monte Carlo simulations which are currently under investi- 
gation [20]. Similar perturbative calculations have been performed in the case 
of isotropic lattice [21,22]. The calculation in this paper is an extension to the 
anisotropic lattice. 

This paper is organized as follows. In section 2, domain wall fermion action 
on anisotropic lattices is given. In section 3, the free domain wall fermion 
propagator on anisotropic lattices is derived. We also study the dispersion 
relation of the free domain wall fermion on anisotropic lattices. It turns out 
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that even in the free case, hopping parameters of the fermion action have to 
be tuned properly, according to the value of the quark mass, so as to have 
the correct continuum limit for the quark. In section 4, the strategies for 
the calculations of fermion self-energy to one- loop are outlined. In section 5, 
numerical results for the one-loop calculation are presented for various bare 
parameters. The renormalization factors for the quark field and the current 
quark mass are given. We also discuss the renormalization of parameter M5 
within tadpole and mean-field approximation. This parameter has to be tuned 
to the right range in order to maintain chiral properties of the fermion. In 
section 6, we will conclude with some remarks and outlook. Dispersion relation 
of the free domain wall fermion is discussed in Appendix A. Some explicit 
formulae for the one-loop calculations are listed in Appendix B. 



2 Lattice Action for Domain Wall Fermions on Anisotropic Lat- 
tices 



The lattice gauge action used in this study is the tadpole improved action on 
anisotropic lattices: [19,7] 
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where Poi and represents the usual temporal and spatial plaquette vari- 
able, respectively. Rij and R^q designates the 2x1 spatial and temporal Wilson 
loops, where, in order to eliminate the spurious states, we have restricted the 
coupling of fields in the temporal direction to be within one lattice spacing. 
The parameter u s , which is usually taken to be the fourth root of the average 
spatial plaquette value in the simulation, implements the tadpole improve- 
ment. Using this gauge action, glueball and hadron spectra have been studied 
within quenched approximation [19,7,8,23,24,25,26]. 

For the fermion action, we adopt Shamir type domain wall fermions [2]. The 
five-dimensional fermion and anti-fermion field is denoted as ij) X)S and ip x ,s, 
respectively. The fermion lattice action we studied is given by: 



Sj — S j + S m , 

Sf = - £(1 + 3k s + Kt- M 5 )xp X!S ip X: 

x,s 
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+ [^x,sPr^x,s+1 + ^x, S Pl^x,s-1 



+ K t^x,s (! + l^)U^(x)ij x+fl:S + (1 - 1^)UI(X - p,)l(). 



S m = -mJ2 [^x,oPl^x,l s -i + ^x,l s -iPriPx,o 



(2) 



We will also use the notation k = n t , Ki = k s . The left and right-handed 
projectors are: Pl/r — (1 T 7s)/2- For definiteness, we take the extension in 
the fifth dimension to be L s and the coordinate of the fifth dimension is la- 
belled such that < s < L s — 1. Thus, our domain wall fermion action is 
characterized by four parameters: five-dimensional mass (wall height) param- 
eter M 5 , temporal hopping parameter n t , spatial hopping parameter k s and 
current quark mass parameter m. 2 . 



3 Free Domain Wall Fermion Propagator on Anisotropic Lattices 

In this section, we will briefly review the results for the free propagator of the 
domain wall fermions on anisotropic lattices. Free domain wall fermion prop- 
agator on an isotropic lattice has been discussed in the literature [2,27,21,22]. 
Our notations follow those in [2,22] and we have made modifications to the 
anisotropic lattice where necessary. We have also tried to keep all relevant 
mass terms (bare quark mass m and the residual mass m r ) in the free domain 
wall propagator. This is useful in practice since we would like to apply the 
anisotropic lattice formalism to massive quarks (charm) and the extension in 
the fifth dimension is usually not very large in practical simulations. 

3. 1 The free propagator 

After performing the four- dimensional Fourier transform, the free domain wall 
fermion matrix appearing in action (2) is given by: 



where the function b(p) and the notation f are defined as: 

2 Note that there is also an additional hidden parameter in the theory, namely the 
extent of the fifth dimension: L s 



D*J(P) = I-Kp) + if ) $ss> + [PrS s +i, s ' + P L S.-i,A 
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b(p) =1-M 5 + J2 Si 1 - C0S Pm) 



(5) 



The matrices VK^p) appearing in Eq. (4) are defined as: 



[^ + (P)]L = W (P) , = -Kp) S ss> - m6 La - ll8f 6 afi + <S._i,« 



(6) 



With this notation, one easily works out the second order operator il^°\p) 
which is defined as: 



>(°)/ 



p 2 + W + W~P R + W~W + P L . 



(7) 



Note that the operators £l^J R (p) = p 2 + W Zf W ± have trivial Dirac structure 
and they are related to each other by: 



^l\p) = nk» 



L a — s— IX,,— s'— 1 



The matrices W~W + and W + W~ are referred to as the mass matrices for the 
right and left-handed fermions, respectively. Free fermion propagator is then 
expressed in terms of the inverse of the second order operators as: 



S F (p) = [D^(p)]\ 1 = [ J D (0) (p)] t ■ [^ (0) (P) 

-if +W + {p) 



= [-if + W (p) 
where we have used the notation: 

Gf /R {p) = [ni%(p) 



P L Gf , 



(9) 
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The explicit results for G R \p) and G^'(p) are found to be: 
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where G^°°^(p) is the corresponding Green's function for infinite fifth di- 

J s,s' 

mension, i.e. when the boundaries are absent. The other coefficients in this 
equation are all functions of the four-momentum p, the boundary mass pa- 
rameter m and the extension of the fifth dimension L s . The explicit formulae 
are: 



G {co \p) 



Be -*c\ss>\ B 
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2b sinh olq 
A± = A- 1 B(1 - m 2 )(e TaG - b) , 

A m = -A _1 m + A~ l Bm r (be~ a ° - 1 + m 2 (l - be aa \ 



A = e 2a ° (b - e- a °) (l + m 2 m 2 r ) + (m 2 + m 2 ) (e a ° - b) 



+ 2mm r b(e 2aa - 1) . (12) 
Here the (momentum-dependent) parameters «g and m r are defined as: 



coshac(p) = 



26(p) 



(13) 



In the perturbative calculation of fermion propagator, it is useful to choose a 
basis such that the mass matrix of the free fermion propagator is diagonal [22] 
in the fifth dimension. In other words, one seeks for unitary matrices and 
V (0) such that: 



(M 2 ) S 6 SS , , IV^W+W-V^} , = (M 2 ) S 5 SS , , (14) 



where (Mq) s for s — 0, 1, • • • , L s — 1 are the corresponding eigenvalues. 3 Note 
that due to the symmetry relation (8), the spectrum of W~W + and W + W~ 
are identical. However, the unitary transformation matrix is not the same 
as V 1 - -*, instead, they are related by: 



vi3 = u {0) 



s,L s —l—s' 



(15) 



The unitary transformation matrices Li^ and are obtained by finding 
all linearly independent eigen-modes corresponding to all possible eigenvalues 
(Mq) s in Eq. (14). The properties of these eigen-modes can be easily obtained 
from corresponding formulae in the isotropic case [21,22]. This set of basis is 
also useful when we discuss the renormalization factors of the chiral mode. 



The suffix in (Mq) s stands for free theory. 
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3.2 Dispersion relation of free domain wall fermions 



By inspecting the low-momentum behavior of the free domain wall fermion 
propagator, one verifies that there exists a left-handed and a right-handed 
chiral pole in the corresponding fermion propagator near the two domain walls 
at s = and s — L s — 1, respectively. These chiral poles reflect the existence 
of the corresponding chiral mode bounded at these walls [2], as can be easily 
verified by solving the corresponding Dirac equation for D^°'(p). The only 
complication that arises in the case of anisotropic lattice is that, in order 
to restore rotational symmetry for small lattice momenta, one has a relation 
among Xi K t and k s when other parameters (e.g. M 5 , m and L s ) are fixed. In 
the massless limit, i.e. both m and m r being zero, this relation simplifies to: 
X = ify/fts- For massive fermions, however, the relation is quite complicated 
and we refer the reader to Appendix A for the explicit formulae (c.f. Eq. (A. 5) 
and Eq. (A. 9)). Here, we will only outline the strategy to derive this relation 
and show some general features of it in figures. 




Fig. 1. The pole mass of the free domain wall quark, thq, measured in l/a s unit as 
a function of the propagator bare mass parameter mp = (1 — b(0) 2 )(m + m r ) for 
three values of the anisotropy parameter x- Three curves corresponds to x = b 3 
and 5, respectively. We have set k s = 1, M5 = 0.5 and L s = 8 in this plot. 

To obtain the dispersion relation for the free domain wall fermion, one searches 
for the pole of the free domain wall fermion propagator. The position of the 
pole in terms of four-momentum is given by: p^ = (iE p , p). This then gives 
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Fig. 2. With the same set of parameters as Fig. 1, the value of Kt/x is shown as a 
function of the pole mass itiq for \ = 1, 3 and 5. It is seen that in the massless 
limit, the expected relation: x = K t/^s is recovered for all x- 



the dispersion relation of the fermion: E p as a function of the three- momentum 
p. The pole mass of the fermion is identified as: rriQ = E p= q = Eq. 

In Fig. 1, we have illustrated the pole mass of the massive quark tuq (measured 
in l/a s unit) as a function of the bare propagator mass parameter mp defined 
as: 



m P = (1 - b(Q) 2 )(m + m r ) 



(16) 



for three values of anisotropy \i namely x = 1> 3 and 5, respectively. Note 
that the residue mass is given by m r = e~ aG ^ mQ ' ' '°^ La , as explained in the 
appendix. The other relevant parameters in this plot are chosen as: k s — 1, 
M 5 = 0.5 and L s = 8. It is seen that the pole mass rriQ basically depends on 
mp linearly for mp values not too large. The nonlinearity only sets in for large 
bare quark mass values. Also, the dependence of ttiq on mp saturates quickly 
for anisotropy larger than unity. For example, the difference between the case 
X = 3 and x = 5 is barely visible. 



The dispersion relation for the free domain wall fermion is quite complicated. 
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However, for small three lattice momenta: a s |p| <C 1, one finds that: 



E P = m Q + -^ + 0(p 4 ), (17) 



where Mk in is the so-called kinetic mass of the quark. To recover normal 
relativistic dispersion relation for small three lattice momenta, one requires: 
rriQ = Mkin which yields a relation among the pole mass and the other param- 
eters in the theory (see Appendix A for the explicit formulae). For the same 
set of bare parameters as in Fig. 1, we have shown the values of Kt/x as a 
function of the pole mass mq (measured in l/a s unit) in Fig. 2. Note that in a 
quenched calculation, the anisotropy parameter x is m fact fixed by the pure 
gauge sector. Therefore, once the pole mass of the quark (or equivalently the 
propagator mass mp) and other bare parameters are fixed, one has to tune 
n t in such a way that the dispersion relation of the fermion resembles the 
continuum form in the small a s |p| region. Basically one utilizes the relation 
shown in Fig 2 to find out the appropriate value of Kt for a give value of mq. 
Of course, for vanishing mq, one recovers the expected relation: x — i^t/i^s- 



4 Fermion Propagator to One-loop 



In this section, we outline the basic strategies to calculate the domain wall 
fermion propagator to one-loop. Numerical results will be provided in the next 
section. To simplify the presentation here, some explicit formulae for the loop 
integrals are provided in Appendix B. 



4-1 Feynman rules 



First of all, one needs the free propagator of the lattice gauge fields [28]. 
This is obtained by expressing the pure-gauge action (1) (with the tadpole 
improvement factor u s set to unity) into the gauge potential A^(x). One also 
needs to fix to a particular gauge. We adopt the following gauge-fixing action: 



gf 



1 1 

2a g x 



E 



X 2 d t {a t A t ) + Y, d 3^sA j ) 



where a g is the gauge-fixing parameter. The gauge field measure action and 
the Fadeev-Popov ghost action do not enter our one-loop calculation. 
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After performing Fourier transformation of the gauge fields, the quadratic part 
of the gauge action then has the standard form in momentum space: 



where 

M 00 = X 



E 



dH 



4(z)iiV(0A,H) 



(19) 
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(20) 

(21) 

(22) 
(23) 



with lattice momentum defined as: Z M = 2 sin(y). The quantities appearing 
in the above equations are given by: 



(24) 



Using these notations, the free gauge field propagator is expressed as: 

F U {1 P A^X)' 



D^l) = M~l 



(I 2 ) 2 I 



fD(ip,q P *,x) J ' 



(25) 



For simplicity, we set a g = 1 in the following calculation. The explicit expres- 
sions for f^ u and fo maybe found in the literature [28]. 

The vertex functions can be obtained from the corresponding expression in 
the isotropic lattice case, which was given in Ref . [22] . After making necessary 
modifications to the case of anisotropic lattice, we have for the one gluon 
emission vertex: 



Vi(k,p,l;a,n) = -igT a n^ 



7m cos \(P» -k^)+i sin ^(p M - k^ 



(26) 



where k, p and I = —{k + p) are the momenta for the anti-quark, quark and 
the gluon field, respectively; a and \x representing the color and Lorentz index 
of the gluon. The anti-quark, quark and two-gluon vertex is given by: 



10 



V 2 (k,pJ 1 ,l 2 ;a,b, t i,v) = ^{T a ,T b } 



5 



flU 1 



(27) 



where k, p, l\ and 1% are the momenta for anti-quark, quark and the two gluon 
fields with the energy-momentum conservation constraint: k + p + l\ + l 2 = 0. 
/i, v, a and b are the Lorentz and color indices for the two gluons. Expressions 
for other vertices will not be needed in our one-loop calculation. 



4-2 Fermion self-energy 



To one-loop order, two Feynman diagrams contribute to fermion self-energy: 
the tadpole diagram and the half-circle diagram. Thus, the 1PI fermion 2-point 
vertex function reads: 



i E Vfr sin Pfl + W + (p)P R + W~(p)P L - S(p) 



(2 C 



where the fermion self-energy S(p) receives contributions from the tadpole 
and the half-circle diagrams, respectively: 



£(p) = S tadpolc (p) + S 



half —circle 



(P)- 



(29) 



Using the Feynman rules given in the previous section, the contribution from 
the tadpole diagram can be written as: 



s tadpolc (p) = ]_ g 2 Cp £ s( gin ^ + cogp ^ 
1 M 



f 



E i 7/*PA«- ? tad,/i( s > S ') + M tad(s, s') 



(30) 



where the explicit expressions for / t ad,/i(s, s') and M ta( j(s, s') are given by 
Eq. (B.l) in appendix B. Note that both Jtad,^(s, s') and M ta d(s,s') are in- 
dependent of the quark mass explicitly. They are diagonal in the fifth dimen- 
sional index and can be computed by numerical integration. 



For the half circle diagram, the contribution can be written as: 
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— 7T ^ 
rr /r\ f ^ + P\ ■ • ^ +P\ 1 

x S F (l) SiS < x ^ 7^ cos(— — + i sm(— — ^ S 



'/■< 



(x 2 (Po-W 2 + E,fe-/. 



„ 2 ) 2 




>(31) 



Note that although £ tadpole (p) j s proportional to the unit matrix in the fifth 
dimensional space, E^* 1 , ircle (p) is non-trivial due to the fermion propagator 
Sf(1)ss' which is given by Eq. (9) in the previous section. To simplify our one- 
loop calculation, we have neglected the effects caused by the finite extension 
of the fifth dimension in the following calculation, i.e. we have set: m r = 
e -a G L s ^ o in the quark propagator Sf(1) S s'- However, in order to study the 
quark mass effects, we have kept the bare quark mass parameter m non-zero. 

The calculation of E s y _circle (p) is somewhat different in the massive and the 
massless cases. If the fermion is massless, it is well-known that E^, 1 , f_circIe (p) 
contains infra-red divergences [28,22]. If the fermion is massive, no infra-red 
divergence shows up in this contribution. In the massive fermion case, since 
there is no infra-red divergence, the expression given in Eq. (31) can be evalu- 
ated directly by numerical integration methods. In the massless case, however, 
the infra-red divergent part (which can be computed analytically) has to be 
subtracted from the self-energy contribution E^, f - circle (p) before it can be eval- 
uated numerically. For the physical quantities in the massive fermion case, we 
denote them as f(m). For the same quantity in the massless case, we denote 
the subtracted expression as: 

/finite = f(m = 0) - /i g , (32) 

where / represents some physical quantity and f\ og contains the logarithmical 
infra-red divergent part to be removed from f(m = 0). Usually, /i og is taken to 
be the corresponding infra-red divergent contribution in the continuum which 
can be computed analytically. After the subtraction of infra-red divergent part, 
/finite is free of infra-red divergences and can be evaluated numerically. 

For the half-circle contribution in the massive case, we obtain: 



S half-circle (m) = _ • £ ^ (jI+( Jn )]^p R + [/" (m)].,,ft) 

- ([M + (m)} s>s ,P R + [M-(m)] s , s ,P L ) . (33) 



In the massless case, after separating the infra-red divergent part, the half- 
circle contribution are obtained as: 
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+ (^finite,^; S ') + 7 log^( S > P L 

- '{M+ nike (s, S ') + M+ s ( S , S '))P R 

+ ((M fi - nitc (s, s') + Mj^{ 8 , s')) P L ] (34) 

where explicit formulae for the quantities appearing in the above formulae can 
be found in appendix B. 

Combining the tadpole and the half-circle contributions, we may express the 
effective action of the domain wall fermion as: 



r( 2 ) = $(- P ) t 



]T z 7M iv [Z+Pr + Z~P L ) + VTP R + W~P L 



V«35) 



where the matrices and W are given by: 



Z^(s, s') = 8 Sj3 > + / tad)/i (s, a') + [J±(m) 

T±(f\\(„ J\ i /\# „/\ , r?\/r±/ 



W^s, s') = W ± (0)(s, s') + M tad (s, s') + [M^m)]^ . (36) 

in the massive case. The quantities [In( m )]s,s' an d [A^ ± (^)]s,s' are replaced by 
the sum of a finite part and an infra-red divergent part in the massless case: 



Zj(s, s') = 5 SjS , + Itad, M (s, s') + I^(S, s') + Ifinite^Cs, s') 



s') = W ± (0)(s, s') + M tad (s, s') + M± g (s, s') + M± ite ( S , s') . (37) 

The explicit formulae for the tadpole contributions J t ad,^(s, s'), M tad (s, s') and 
the half-circle contributions [I^(m)] SjS >, [M ± {m)] SyS , M^ g (s,s'), M^ nitD (s, s'), 
/ io glM ( s ' J finite, M (X s ') can be found in appendix B. 



4-3 Wave-function and mass renormalization of the chiral mode 



We are concerned with the renormalization of the chiral modes which are 
bound to the walls. To obtain this information, it is better to use a new basis, 
ip d (p), which diagonalizes the mass matrices W in Eq. (36). This new basis 
is related to the old one via: 



i> d M = u s ^p r Mp) + v^Vvb), (38) 
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where the two unitary matrices U and V satisfy 



M 2 S 5 S , S , , 



VW + W V f 



(39) 



Under this new basis, the effective action of the fermion becomes: 



T^=ij d (-p) f 



+ VW + U^P R + UW V ] P L 



(40) 



The unitary matrices W, V and the corresponding eigenvalues (M 2 ) s can be 
calculated to 1-loop level: 



U = (1 + g*UM)U® , V = (1 + 2 V«)V< o > , 
(M 2 ) S = (M 2 ) S + /(M 1 2 ) S . 

where tree-level matrices Wq and Vo are defined in Eq. (14). 



(41) 
(42) 



After diagonalization of the mass matrix, the effective action for the chiral 
mode field ip d (p)L s -i = Xo(p) becomes: 



Xo(-p) 
where 



J2 il»P» (Z + ,„Pr + Z-,»Pl) + W + P R + W-P L 



Xo(p) , (43) 



Z ± ,, = l + g 2 (A 



,tad,/i + I±,v,{ m , 



L.-1,L.-X 



W± = —m P 



Mt(m 



L„-1,L S -1 



(44) 
(45) 



The above formulae work for the massive case. In the massless case, one has 
to modify the right-hand side of Eq. (44) and Eq. (45) accordingly as specified 
by Eq. (32). Note that chiral symmetry ensures that the effective quark mass 
term is renormalized by a multiplicative factor, i.e. there is no additive renor- 
malization. The matrix elements of various I d and M d in the above expressions 
are given explicitly by: 



2jd W (0)j+ W (0) t 2jd = V (0)j- V (0)t 
,J ~r)P P 7 — )P P 
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n 2jd _ 7/(0)1-+ 7/(0)f 2jd _ V (0)t- y(0)f 

a 2rd -u(°)t+ u m o 2 i d -v^r v (0)t 

y J +,finite,At — M J finitc,p M > i) J finite,^ — K J finite, p v 

9 2 m P M d +Md = V^M^U™ , g 2 m P M d _ Md = U^M tad V^ , 
g 2 m P M d + = V (0) M + W (0)t , y 2 m P M^ = U^M~V m , 



g 2 m P M d +Mitc = V^Ml itc U^ , g 2 m P M d _ Mitc 
It is also easy to show that: 



(46) 



(47) 
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-1 


V^— , finite, ii) L a — 1,L S - 


= fd 

"1 finite, £t 


« tad ) Lfl 




-1 


= (Mt itad ) La -i,L s -i 


= <d , 






-1 


= {Mi Wi = 




«iogk 




-1 


= (^-,log)i s -l,i s -l 


= < g , 


(^+,finite)L s 




-1 


= .finite) L s -l,i s - 


1 = ^finite 



(48) 



Thus one obtains: Z A 



■ii 



Z M and 



W- = W. Using the above 



notations, we get for the wave-function and mass renormalization factors: 



- 1 + y 2 [IU* 



W± = —m P 



g 2 (M t d ad + M d (m] 



(49) 
(50) 



in the massive case and 



1 



9 2 (4 



tad,^t 



jd 

log,At 



W± = —mp 



g 2 [M i 



rd 
'tad 



finite, /j, 



< + ^finite 



(51) 

(52) 



in the massless case. 

The tadpole contributions I d ad41 and M d ad can be obtained numerically by us- 
ing the definition (48) and the explicit expressions of itad,/i( s ; s') and M ta d(s, s') 
given in appendix B (c.f. Eq. (B.l)). Using the definition of hi^ and V^°\ we 
get: 



rd 

tad, /i 



-c F - 



2 1 w ix% + Erf 



2\2 v t 1 



fix 



/ 



Id 



(53) 
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and if we define M^ ad as 



-C F - 



E 



2_i (2tt)^ (x 2 l 2 o + Ejl 



'2\2 v ^ 
J 



(Ex + 4-) 



/J 



(54) 



we arrive at: 



Mi 



tad 



2e-7((l-& 2 (0))(l + m 2 e- 2a ))M t ° E 



tad > 



(55) 



where the parameter a is exponential decaying rate for the chiral mode in the 
fifth dimension. Using these equations, numerical values of 7 t d ad and Mf ad will 
be presented in the next section. 

We now proceed to evaluate the half-circle contributions: 1^ and M d . If we 
use the shorthand notation: 



and a similarly notation associated with , we have: 



(G L ) U = B 



{Gr)u = B 



sinh ac — sinh a 
cosh oiq — cosh a 
sinh (Xq — sinh a 



cosh «c — cosh a X> 
(W+G £ ) M = (W^G^y = (e~ a - b)GR + R + R 1 
(W+G L ) V = (W Q -G R )u = (e" Q - b)GR + R + R 2 



(Gr)v = Gl 
(Gl)v = Gr 



where: 



(1 -m 2 



A/"i = (l -m 2 )(e" G -6) + 
N 2 = (l-m 2 )(e~ aG -b) - 
V = (1 + m 2 e" 2a ) [{e~ aG - b) + m 2 (fe~ 2aG - e~ aG )] (e aa - e~ a ) 2 



2 )(e~ aG -b)m 2 e- 2a 

a 



■ m 



! )(e° G 



1 - e- 2 « 
-n-o - ^ ; 7T- 77, 



+ 46m 2 sinh age Q 
b)m 2 e~ 2a + 46m 2 sinh a^e 

' e OG _ 



— a 



1 + m 2 e~ 2a (1 - e~( Q 
- i m me aG (e" aG + 2e 



-Q 



<g))2 



1 - e~ 2a 



A " e 

-/L m m— ^-3^ 
1 + m^e 



i? 2 = -A, 



. m m 



e -(a+a G ) 
1 

1 + m 2 e~ 2a 1 - e-( a+a e) 



'm 2 (l - e - 2 ° G )e- 3Q 

Q - 2e" aG e- 2Q - e~ 3a ) 
m 2 e~ 2a 



*)} 



e -2a 1 

1 - e~ 2a 



1(3 



Using the above notations, after some algebra, l£ is found to be: 



I«{m) = C F 



dH 



+ — sin 2 1^{G l + G R ] 



« At (cosZ M ((e a -b)G R + R 



p , no 
J {x 2 il + Ei 



+ 



S ((e- a - + i? + ^li-^) - 2*; 2 cos 2 l,/2G L 



4(X 2 ) 5 



^0 



x 2 /g + E,/l (x/ ' + MO 



2X 



tf(0 , /f(0/ w (0 



+ 



- 2*kJ (sin 2 / At /2G fl + cos 2 IJ2G L 
+ X> 2 ( sin2 ^/2^R + cos2 



4(X 2 ) 5 



x 



fid 



sin 2 Z M 

(x% 2 + Ei §) s 



J /d(0 



"/MO 



(56) 



In the massless case, what needs to be calculated numerically is the subtracted 
finite part ifi nite „ defined as: 



(1-6 2 (0))(1 



-2a\ 



X 



(1 -b(0)e~ a ) 2 

dH 4x(x a )'' rf, ?^-(x 2 ig + E i ?)) 



(2tt) 



(x 2 /g+E ? a 



(57) 



The infra-red divergent part that is subtracted from 7^(0) is denoted as 
and can be obtained analytically by using Eq. (B.8) and the definition (48). 

For the contribution M d , we have: 



(^%G L U^) L ^ 1>L ^ 
(V^GnU^^r, 



{U^G R V^) Ls ^ x 

{U^G R V^) Ls ^r X 



(5* 



One also verifies that: 



qM _ QM 



-mpG 



M 
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—rap 



(l-e- 2a )B 
(1 -6 2 (0))(l + m 2 e" 2a 
A- 1 



(1 _ g-(a G +a))2 

We also find: 



2e~ a (l - m 2 ) (cosh a G -b) 



]_ _|_ e -( a G+&) g-Q 

1 1 _ e -("G+«) i _ e _2a 

1 + m 2 e~ 2aS 



B 



(59) 



(V^W + G L U^) Ls ^ Ls ^ = (wwrc^.^ 

= -m P ((e" Q - + < f + i?f ' 

= -m P ((e" Q - b)G M + itf + R 2 



= —mp 
x | 1 + 



1-e 



-2a 



(1 -6 2 (0))(l + m 2 e- 2Q ) 
2e- a (e~ a -6(0))\ 



1 - e 



-2a 



where 



M 



(1 - e-^BA- 1 



(1 - 6 2 (0))(1 + m 2 e" 2Q ) (1 - e-( Q G+«))2 

x (1 - m 2 ) [(e~ QG - e~ a )(e aG - b) + (e QG - e- Q )(e- aG - 6) 



1-e 



-2a 



(1 -6 2 (0))(l + m 2 e- 2a ) 

f l + 2e-( QG+Q ) 1 / (1 -m 2 )(e aa -b) 
x { B- ; — -T- + — LB- 



1 _ e -(a G +a) £ 



1 _ e -( a G+«) 



-CtQ 



+ 



e~ a )m 2 e~ 2a + e aG - e~ a 



(1 _ e -(a G +a)); 



m 2 e a \ 
1 — e -("o+«) / 



-2a 



pM _ 1 e 

""■2 - ^ !.2^\W1 , „,2„-2^ 



(1 - 6 2 (0))(1 + m 2 e 
x IB 



1 _ e -("G+«) 



A 



+ 



i"e~ 2a + 2e~ {aa+a) t 1 f (1 -m 2 ){e aa - b) 

1 — e -(«Q+a) 

m 2 e~ a \ 

(a G +a) J 



m 2 e~ 2a 



(e aG - e~ a )m 2 e~ 2a + e~ aa - e 



(1 _ e -(a G +a))2 l_ e 

Using the above expressions, the quantity M d is found to be: 



M d = g 2 C F / E "J f cos2 W2 ((^° - &)G M + < + < 
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- sin 2 IJ2 ((e" a - b)G M + + R% ! ) + k m sin 2 l^G M 
x + Ml ) . (60) 
Similarly, in the massless case, the expression for Mg nitc reads: 



M finite = M (m = 0) - 2, C F J — ( X 2 f g + S . f 2 )2 (61) 

Using the definition (48), the subtracted part, i.e. Mf Qg , is given by Eq. (B.ll) 
in appendix B. 

Once the values of J t d ad , M t d ad I*(m) and M d (m) (4^^ and M^ nitc in the 
massless case) are obtained numerically, Eq. (49) and Eq. (50) (Eq. (51) and 
Eq. (52) in the massless case) then give all the relevant renormalization factors 
for the chiral fermion mode to one-loop order. 



5 Numerical results for the one-loop calculation 



In this section, we present numerical results of our calculation of the fermion 
propagator to one-loop. We will present the results for the wave function 
renormalization for the light fermions, the mass renormalization for the mas- 
sive quark and the shift in the parameter M 5 which is important in order to 
have chiral modes. 



5. 1 Wave function renormalization for the domain wall fermion 



In this section, we concentrate on the wave function renormalization constant 
associated with the chiral mode of the domain wall fermion. We will be dis- 
cussing the massless case and the massive case separately. 

For the massless case, the formula for the wave-function renormalization con- 
stant is (c.f. Eq. (51)): 

Zfj, = I + g 2 (Ald,^ + I\og,fi + -^finitc,^) J (62) 



where J^j comes from tadpole diagram and I\ og ^ while /fi nite /J designates the 
logarithmic divergent and finite part of the half circle diagram, respectively. 
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Table 1 

Values of 7^ ad for three different values of \ are listed. Due to symmetry, I\ 
lL, a fori = 1,2, 3 and 7^0 = 7^. 



x = i 


X = 3 


X = 5 


-'tad,! 


tad,s 


J tad,t 


Td 

tad,s 


-'tad,! 


Td 

tad,s 


-0.090487 


-0.088871 


-0.014293 


-0.122173 


-0.005399 


-0.127386 



The quantity If og ^ as defined in Eq.(48) and Eq. (B.8), can be evaluated 
analytically. The value of Ii nite ^ is obtained from Eq. (57) by numerical inte- 
gration. Similarly, using Eq. (56) one finds the value of I^{m) in the massive 
case once the bare parameters in the fermion action are given. The value of 
Itadfi i s obtained by performing the integration appearing in Eq. (B.l) numer- 
ically. 

In Table 1, we have listed the tadpole diagram contributions to / t d adM for 
several values of anisotropy. Due to symmetry, we have denoted 7 tad i = J t d ad s 
for i = 1,2,3 and J^ d = J tad t . The values in Table 1 are obtained for \ — 1,3 
and 5. Note that the tadpole contribution / tad;At only depends on the anisotropy 
parameter \ an d is independent of other parameters in the fermion action. 
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-0.001 - 
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finite ti-0.01 0- 
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-0.015- 
-0.016- 
-0.017- 
-0.018- 
-0.019- 



• 





A * * 



4 - 
▼ ▼ 



▼ ▼ 



▼ ▼ 



(X=1) 

(x=1) 
(X=3) 
(X=3) 
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0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 3. The values of I^ nitc „ as a function of M§ for three values of anisotropy 
X = 1, 3, 5 at k s = 1.0, K t = xk s . 

For the massless case, the finite part of the half-circle diagram contribution, 
i.e. /|aite,u) are calculated using numerical integration when the parameters of 
the theory are fixed. In Fig. 3 and Fig. 4 we have plotted the values of Ifl n i te ,u as 
a function of the parameter M 5 for three values of the anisotropy (x = 1, 3, 5) 
and two values of k s (k s = 1.0,0.5). Since /g niteM always corresponds to the 
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finite 
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Fig. 4. Same as in Fig. 3 but with k s = 0.5. 

massless case, we have set n t = x K s in these plots. The infra-red divergent 
contribution from the half-circle diagram can be obtained from Eq. (B.8) once 
a infra-red scale is given. 

For the massive case, one needs the values of I^{m) which is obtained by 
numerically integrating Eq. (56). The values of I^{m) are shown in Fig. 5 
through Fig. 8. Similar to the massless case, the corresponding values of I^(fn) 
are plotted as functions of M 5 for three values of anisotropy: x = 1)3,5 and 
two values of k, s (k s = 1.0,0.5). 



5.2 Mass renormalization for the domain wall fermion 



For the mass renormalization factor given by Eq. (52), we obtain: 
Z m = 1 + g 2 (Mt d + Mf OR + M| nite ) , 



(63) 



for the massless case and for the massive case we have: 
Z m = l + g\Mi d + M d ). 



(64) 



For the tadpole contribution M t ^ d , using notation (55) and setting k s = 1.0, 
we have evaluated the tadpole contribution M t ° ad for three values of anisotropy 
and the results are listed in Table 2. The values of M t ° ad listed in Table 2 are 
obtained at: n s = 1, M 5 = 0.5 with the values of K t being tuned accordingly. 
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Fig. 5. The values of I^{m) as a function of M5 for three values of anisotropy: 
X = 1, 3, 5 at current quark mass m = 0.3. The parameter k s = 1.0 and the values 
of Kt are adjusted accordingly, as discussed in Sec. 3. 



Table 2 
Values 
to those at k s 



Values of A^ t ° ad for % = 1, 3, 5 and m = 0, 0.3, 0.5. The values listed here correspond 



1, Mr = 0.5 and the va 



ues of Kt are tuned accordingly. 





x = 1 


X = 3 


X = 5 


m = 0.0 


-0.357099 


-0.409400 


-0.409156 


m = 0.3 


-0.355505 


-0.408316 


-0.408381 


m = 0.5 


-0.352263 


-0.406204 


-0.406878 



For the half-circle contributions in the massless case, values of Mg nite are 
presented in Fig. 9 and Fig. 10 where we have shown the values of Mg nite 
for k s = 1.0 and k s = 0.5, respectively. In each figure, three values of the 
anisotropy, namely x — 1> 3, 5, are taken to evaluate Mg nite as a function of 
M 5 which is taken to vary between zero and unity. 

For the massive case, we evaluate M d (m) directly by numerical integration. 
In Fig. 11 and Fig. 12, we have shown the values of M d (m) as a function of 
M5 at m = 0.3 and m = 0.5, respectively. We have also taken three values for 
the anisotropy parameter in each case. 
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Fig. 6. Same as in Fig. 5 with k s = 0.5 and m = 0.3. 



5.3 Mass shift for the parameter M5 in mean-field theory 



At the tree level, it is known that the domain wall fermion preserves chi- 
ral properties only when the parameter M5 is tuned to the right interval ( 
< M5 < 1 for free domain wall fermion). Since this parameter enters the ac- 
tion via a five- dimensional mass term, one expects that it receives substantial 
additive renormalization. At one-loop level, it is important to check that, after 
renormalization effects are taken into account, whether this so-called "chiral 
window" for M 5 remains there or not. In the case of isotropic lattices, it has 
been shown that the window remains in perturbation theory and one expects 
chiral properties of the fermion as long as M 5 lies in the right interval. We 
wish to check this point for the anisotropic lattices. 

In perturbation theory, as expected, the main contribution for the shift of 
the parameter M5 comes from the tadpole diagram. To simplify the analysis, 
we will only discuss the shift of M 5 due to tadpole diagrams. Within this 
approximation, one finds that the parameter M5 is replaced by: 

M 5 = M 5 - Ml - U t ) + 3k,(1 - U s )) = M 5 - 5(M 5 ) , (65) 



where the tadpole parameters U t and U s are given by: 
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Fig. 7. Same as Fig. 5 with current quark mass parameter m = 0.5 and k s = 1.0. 

1 J dH 1 - / 00 

here g 2 and x are the "boosted" coupling and the "boosted" anisotropy [29] 
which are defined as: 



9 2 = ^U!U t , X = Xo^. (67) 

Eq. (65) suggests that in perturbation the chiral window remains but it is 
shifted. 

We will only show results at zero current quark mass m = in which case 
we have: KtJ k s = x- To get a feeling how large the shift in the parameter 
M 5 can get, we take a particular value of (3q. At first, the values of U% and 
U s are set as U t — U s — 1. By making use of the formula (67), the values of 
g 2 and x are obtained and can be substituted into formula (66) to obtain a 
new set of values for U t and U s . The above procedure is repeated iteratively 
until the values of U t and U s become stable. The ultimate values for U t and U s 
in perturbation theory and the corresponding shift in M5 thus obtained are 
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Fig. 8. Same as Fig. 5 with current quark mass parameter m = 0.5 and n s = 0.5. 

listed in Table 3 for (3q = 3.0 at five different values of xo- 
Table 3 

The values of tadpole improvement parameters Ut, U s and corresponding values for 
5M$ as evaluated from Eq. (65) at xo = 1 to xo = 5. The bare gauge coupling is 
Po = 3.0 and the parameter k s = 1. 





Xo = 1 


Xo = 2 


Xo = 3 


Xo = 4 


Xo = 5 




0.889580 


0.956099 


0.977525 


0.986559 


0.991122 




0.891130 


0.861399 


0.851613 


0.847194 


0.844862 


S(M S ) 


0.437222 


0.494908 


0.503901 


0.504587 


0.503253 



Another way to investigate the renormalization of parameter M5 is to use 
mean-field approximation. In this case, the domain wall fermion action is 
rewritten as: 

[<SW] S , S ' = (~Kp) + * 2 K yPifln sinp M ^J 8 StS > 

+ P R 8 a+lt8 > + P L 5 s -i,s' , (68) 

where are pure numbers 4 and b{p) is defined as: 
4 We use the notation: Ut = Uq and U s = U{ for i = 1, 2, 3. 
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Fig. 9. The values of M^ nitc as a function of M5 are shown for % = 1, 3 and 5. Other 
parameters are: k s = X K s, m = 0. 
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Fig. 10. Same as Fig. 9 except that k s = 0.5. 



b(p) = 1 - M 5 + ^2 ac m (1 - Up cospp) 



(69) 
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Fig. 11. The value of M d as a function of M5 is shown for \ = 1, 3 and 5. Other 
parameters are: k s = 1, m = 0.3. 

Following similar calculations as for the free domain wall propagator, one can 
define an effective mass for the fermion in the mean-field approximation. We 
will denote this mass parameter as mp. 

In Fig. 13, we show the effective mass mp of the chiral mode as a function 
of the parameter M5 for five values of xo- F° r simplicity, we have set k s — 1 
in this figure and the extent in the fifth dimension is taken to be L s = 10. 
It is seen from the figure that, within a certain interval, the chiral window 
remains for all Xo- This has been checked in the case of isotropic lattice which 
roughly corresponds to the red curve in our figure. We see from the figure 
that this qualitative feature remains valid in the anisotropic lattice case for 
almost all values of anisotropy Xo- I n fact the shape of the curves is quite 
insensitive to the value of Xo-, as is seen from the figure. The size of the chiral 
window remains more or less unchanged but the starting and ending position 
is shifted. The results discussed above suggest that, as long as the parameter 
M5 is tuned to the right interval, massless chiral modes are stable against 
perturbative quantum fluctuations. This result also serves as a guidance for 
the choice of the parameters in our future numerical simulations. 

We have also investigated the dependence of the chiral window on the size of 
the fifth dimension L s . In Fig. 14, we have plotted the effective mass of the 
chiral mode as a function of M 5 at two different values of L s . The anisotropy 
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Fig. 12. The value of M d as a function of M5 is shown for \ = 1, 3 and 5. Other 
parameters are: k s = 1, m = 0.5. 

parameter is set to Xo = 5 while all other parameters are kept the same as 
those in Fig. 13. It is seen from the figure that the chiral window for the smaller 
lattice (L s = 8) becomes narrower than that of the larger lattice (L s = 16), 
as expected. 



6 Conclusions 



In this paper we have studied the domain wall fermions on anisotropic four- 
dimensional lattices. We have analyzed the free domain wall fermions, showing 
that the hopping parameters have to be tuned according to the current quark 
mass parameter in order to restore the correct dispersion relation for the chi- 
ral mode of the domain wall fermion. Using lattice perturbation theory, the 
domain wall fermion self-energy is calculated to one-loop order. We obtained 
the wave-function renormalization constant and the mass renormalization 
constant Z m for the chiral mode. It is also verified that the chiral property of 
the fermion remains when quantum corrections are considered. We have also 
estimated the range for the parameter M 5 in order to maintain a chiral mode. 

From our study, it seems that domain wall fermion on anisotropic lattices is 
applicable in practical Monte Carlo simulations. The only parameter that 
has to be tuned is the hopping parameter n t when other parameters are 
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0.12 n 




0.5 1.0 1.5 2.0 2.5 
M s 

Fig. 13. The effective mass of the chiral mode, mp, as a function of the parameter 
M5 at m = 0.0, L s = 10. Five curves with different colors correspond to different 
values of xo- The values of Ut and U s are taken to be those in Table 3 at (3q = 3.0. 



given. The tuning of this parameter can be done either perturbatively or 
non-perturbatively. The non-perturbative tuning requires the measurement 
of physical hadron dispersion relations or quark dispersion relations by Monte 
Carlo simulations. An exploratory numerical study is now under way and we 
hope to come up with some results in the future [20]. Compared with similar 
situations for the Wilson fermions, anisotropic lattice domain wall fermion 
appears to be less complicated since there is no need to add clover terms (or 
other dimension five operators) due to better chiral properties. Of course, this 
conclusion relies on the fact that the residual mass effects due to the finite 
extent of the fifth dimension to be small. We have tried to estimate this ef- 
fect with mean-field approximation but this issue definitely should be checked 
further in future numerical studies. 

It should be straightforward to generalize the results we obtained in this paper 
to the overlap fermions on anisotropic lattices. Another topic that is worth 
studying is the renormalization factors for quark bilinear operators. The study 
of these problems are now under investigation. To summarize, we anticipate 
the anisotropic domain wall fermions to be helpful in the study of the heavy- 
light hadronic systems, exotic hadrons (hybrid) with light quarks and light 
hadron-hadron scattering processes and we hope to come up with some ex- 
ploratory numerical results in the future [20]. 
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Fig. 14. The effective mass rap of the chiral mode as a function of M5 for two 
different values of L s . The red and black curve corresponds to L s = 16 and L s = 8, 
respectively. Other parameters are the same as those in Fig. 13. 

A Dispersion relation of free domain wall fermions 



In this appendix, we briefly outline the derivation of the dispersion relation for 
the free domain wall fermion on anisotropic lattices. The dispersion relation 
is such that the four- momentum p = (iE p , p) is a zero of the function: 

A = (1 + mm r ) 2 be 2aa - (1 - m 2 )(l - m 2 r )e aG - (m + m r ) 2 b . (A.l) 



In principle, the parameter m r = e~ a ° Ls also depends on the four-momentum. 
For the moment, we assume that the extension in the fifth dimension is large 
enough such that we can ignore the momentum dependence of m r in the 
study of the dispersion relation. This amounts to taking m r = e~ aG ^ Ls . If we 
denotes: 

mi = (1 + mm r ) 2 , m 2 = (1 — m 2 ) (1 — m 2 ) , m 3 = (m + m r ) 2 . (A. 2) 



Then, the equation for the pole yields: 
- p 2 = Mi [(1 + b 2 + p 2 ) 2 - Ab 2] 



(A.3) 
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where M\ = mim^/ml. The complete dispersion relation turns out to be quite 
complicated. But if the lattice three-momentum |p| <C 1, we obtain: 



E r 



E + 



2M 



kin 



+ 0(p 4 



(A.4) 



where E will be identified as the pole mass of the quark: E = rriQ and M^ in 
is the so-called kinetic mass. After some calculations, the equation satisfied 
by Eq = iTiQ is found to be: 5 



A cosh 2 E + B cosh E + C = 



(A.5) 



where 



A = n 2 



(6(0) + «t) 



Mi 



5 = 4(6(0) +«t)«t (6(0) + «t 



ft, 



Mi 



(6(0) + ^) 4 + 2(6(0) + ^) 2 K 2 



(«? + I) 2 



We can also find out the kinetic mass term with the result: 



2M, 



kin 



dE 
dp 2 



p=0 



Ai cosh 2 E + Bi cosh E + C x 
(2A cosh E + B)siuhE 



Mi (A.6) 



(A.7) 



where 



A 1 = 4(6(0) + n t )K 2 t K s M 1 



By 



3(6(0) + K t f + 2(6(0) + K t )K s + K 2 t - 1 K t K s M x 

Ci = k 2 + 2k s [(6(0) + re*) ((6(0) + ^) 2 + K 2 - l) 
+ k s ((6(0) + K t ) 2 + «: 2 + l)] Mi 



(A.8) 



In order to have the usual energy-momentum dispersion relation, one imposes 
the condition: Mk in = mq which yields the relationship between \ an d i^ti K s 
as (restoring the lattice spacing explicitly): 



X 



(2 A cosh E + B)sinhE 
2(Ai cosh 2 E Q + B x cosh E Q + d)E 



(A.9) 



5 Note that the quantity Eq appearing in this equation is the dimensionless quan- 
tity: atEo = a t mQ = (a s mq)/x, if we restore the dimension for Eq. 
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where Eq = tuq is obtained by solving Eq. (A. 5). Eq. (A. 5) and Eq. (A. 9) thus 
establish the relation among the pole mass tuq, anisotropy \ and the other 
bare parameters of the theory for free domain wall fermions on anisotropic 
lattices. 

In the above discussions, we have neglected the momentum dependence of 
m r = e~ aG ^ p ' Ls . If we keep this momentum dependence, the pole position 
can be solved numerically. For any given set of parameters, we start with the 
choice m r = e~ aG ^ La and solve for niq in Eq. (A. 5). Then the obtained value 
of Eq = rriQ is substituted into m r = e - a a(^ Q ,o,o,o)L s an( j a new va j ue Q f £j Q 

is thus obtained by solving Eq. (A. 5) again. This procedure can be iterated 
until a stable set of values for rriQ and m r = e~ aG ^ mQ ' ' '°^ L3 is obtained. The 
figures in the main text are obtained using this method for a given L s = 8. 



B Some Explicit formulae for loop integrals 



In this appendix, we list the explicit expressions for various loop integrals 
which enters the fermion self-energy discussed in section 4.2. 

First, the tadpole contributions itad,,u and M ta d are given by: 



1 f rl 4 l k * 

M tad ( Sj s') = -g 2 C F - / f4V , " x (P x + l—)5 ss ,. (B.l) 

Note that these tadpole contributions are diagonal in the fifth dimension. It 
is also noted that /tad,/j and M ta d are independent of the current quark mass 
parameter m explicitly. 

Now comes the formulae for the half-circle diagrams. In massive case, we have: 



[It(m)] s , s , = g 2 C F J 



dH 



(oobI^W-Gr + W + G L )(s, s') 



+ K M sin 2 lfj,(G L + G R )(s, s'] 



Xl 



2 , /""(Q 
^ fo(l) 



+ 



(x 2 l 2 + Ej W 

(^(W'Gr + W + G L ) - 2kI cos 2 IJ2G^\ (a, s') x 
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MO 



2x 



4(x 2 )^° 

2 X kI (sin 2 Z M /2G ii + cos 2 lj2G L ) (s, s' 



fn(l) fh(l) 



K 2 V (sin 2 l„/2G R + cos 2 l„/2G L ) (s, s' 



4(x 2 )^° _ {x p , no, /;(0 , f?(i)r(i) 



X^ + E^ 2 



f D (l) ' fnil) 



fid) 



sin I 



(x 2 l 2 o + Ej Z 2 ) 2 



[4( m )]«,«'=0 2c i 



d 4 z 



+ k m sin 2 l^(G L + Gh)(s, s') 



x/ 2 + 



/ MM (0 

foil) 



(x 2 l 2 + Ej Z 2 ) 2 



+ 



-K M (WGf fl + - 2«J cos 2 Z m /2G r ) (s, a') 



(B.2) 



2\5 



4(X 2 ) 



fi() 



2X 



MO /SCO 



- 2*«J (sin 2 IJ2G L + cos 2 Z M /2G fi ) (s, s') 
+ 5> 2 (sin 2 Z„/2G L + cos 2 l v /2G R ) (s, s' 



x 



X 



x 2 Z 2 + E,Z 2 

sin 2 Z^ 
(x 2 § + Ej i 2 ) 2 



foil) ' fn{l) 



fUi) 



[M + (m)} s , s , = g 2 C f 



•3 3 ' 

dH 



E«2 



(2tt) 4 z 



cos 2 lj2{W + G L ){s,s') 



(B-3) 



- sin 2 L/2(W-G R ){s, s') + - Klt sin 2 Z M (G L + G R )(s, s') 



x 



(x 2 ^ + E, i 2 ) 2 



(B.4) 



[M-(m)] v = ^Ci 



'3 3 ' 

dH 



(2tt^ 4 



cos 2 Z„/ 2 (^~ G *)( s > s ') 



- sin 2 Z„/2(W+G L )(s, s') + -«„ sin 2 Z M (G L + G fl )(s, s') 
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X 



(x 2 i 2 + Ej I?) 2 



(B.5) 



where we have used the following notations: 



sin l^dpn 



df D (p-l) 



p=0 



sin l^dpn 



(B.6) 



p=0 



It's easy to verify that M ± (m) is proportional to mp = (1 — b 2 (0))m. Thus, 
we may define: 



[M ± {m)} s ^ = (1 - b 2 (0))m{M±(m)] f 
In massless case: 



O s > s ') = Y^9 2 C F (C ± ) s , sl (hi(7T 2 ) + 1 - \n(p 2 ) + In k 2 s 

iLt e A s ^ s ') = K( m = °)jV - g 2 c F (c + ) s , s , 

f dH 4 X ( X 2 )<W^(7r 2 - (X% 2 + gjffl) 

X J (27T) 4 (X^O + Ej ?) S 

f dH 4 X ( X 2 )^/^ (7r 2_ (x2/ 2 + £ , / 2 )) 

X 7 (2vr)4 ( X 2 / 2 + Ej Z 2 ) 3 

1 



M io g (^ = -^9 2 C F (CM) S , S , (ln(7r 2 ) + 1 - ln(p 2 ) + In* 2 
MLte(s, s>) = m P {[Af +(m = 0)] v 



fimtc( S > S ') = m p{[ M (™ = 



' /• ci 4 / 4 X fl( 7 r 2 -( X 2 / 2 + E J / 2 )) ] 
J (2ny (X% 2 + E,/ 2 ) 2 J 



-9 2 C F 



mp 



where (C±) S)S ' are given by: 



:< /• gj 4 X g(7r 2 -( X 2 / 2 + E^ 2 )) l 

"7(2^ (x% 2 + E^ 2 ) 2 J 



(B.7) 



(B.8) 

(B.9) 

(B.10) 
(B.H) 

(B.12) 
(B.13) 



(Cf) s , s , = (l-6 2 (0))m 

+ 6(0) (£W + {ti) s , t e~* {L '- x ~^\ 
(C M ) SyS , = (1 - 6 2 (0))m koe-" ^" 1 -^ 



(B.14) 
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+ 6(0) E W-(0) s , t e- ao{L °- 1+t - s,) 



(B.15) 
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